Measuring photon anti-bunching from continuous variable sideband squeezing 
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We present a technique for measuring the second-order coherence function g*-^^ (r) of hght using a Hanbury- 
Brown Twiss intensity interferometer modified for homodyne detection. The experiment was performed entirely 
in the continuous variable regime at the sideband frequency of a bright carrier field. We used the setup to 
characterize g*^^^ (r) for thermal and coherent states, and investigated its immunity to optical loss. We measured 
(;'^'(t) of a displaced squeezed state, and found abest anti-bunching statistic of (;'^'^'(0) = 0.11 ± 0.18. 

PACS numbers: 03.65.Ta 42.50.Xa 42.50.Ar 42.50.Dv 



Fifty years ago, Hanbury-Brown and Twiss (HBT) first 
demonstrated an optical intensity interferometer Since 
then, HBT interferometry has been applied to diverse ar- 
eas such as condensed matter physics, atomic physics, and 
quantum optics ; and has become a powerful mea- 

surement technique in astronomy, and high-energy particle 
physics jsll^. From a historical perspective, HBT reported 
coiTelations in the intensity measured at two locations, from 
light emitted by a thermal source. The effect was interpreted 
as being either a manifestation of classical wave theory, or 
due to bunching in the arrival time of photons. Such corre- 
lations were generalized to «th-order by Glauber in a com- 
prehensive quantum theory of optical coherence |7], with the 
second-order coherence g^'^'^ corresponding to the measure- 
ment made with a HBT interferometer Curiously, the the- 
ory predicted that certain states of light would exhibit a pho- 
ton anti-bunching effect, which is the tendency for photons 
to arrive apart from one another This is a non-classical phe- 
nomenon which violates the Schwartz inequality |4|. Photon 
anti-bunching has been observed in resonance fluorescence, 
conditioned measurements of parametrically down-converted 
light, pulsed parametric amplification, quantum dots, and 
trapped single atoms/molecules \^ % J^O, 121 . Recent 
experiments have probed the spatial and temporal second- 
order coherence functions of atomic species in BEC, and atom 
lasers 1 13, 14|. All of these experiments have relied upon the 
ability to detect individual particles in a time-resolved mea- 
surement. 

In this Letter we apply a technique for measuring the 
second-order coherence of optical fields, that complements 
previous studies and provides a link between discrete-variable 
(DV) and continuous-variable (CV) quantum optics. Our 
scheme is based on the HBT interferometer, but uses homo- 
dyne detection in each arm, to make CV measurements of the 
quadrature amplitudes over a range of sideband frequencies. 
The second-order coherence 5'^' is then constructed from 
the set of four permutations of the time-averaged correlations 
between the amplitude/phase quadratures. At no point is it 
necessary to make time-resolved detections of single-photons 
fl5]. Homodyne detection offers the advantage of high band- 
width, and excellent immunity to extraneous optical modes. 
We used the scheme to measure the temporal second-order 
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FIG. 1: (Color online) Schematic of the experimental setup. OPA 
optical parametric amplifier, A/2 half-wave plate, PBS polarizing 
beam-splitter, x:y beamsplitter with transmission x, H1/H2 homo- 
dyne detectors, AM amplitude modulator, lo function generator, 
WNG white noise generator, 8 mixer, LP low pass filter. 



coherence function g'^^ (t) of a displaced squeezed state. 

In contrast to most CV experiments involving squeezed 
light, is the realization that weaker squeezed states can ex- 
hibit a greater anti-bunching effect. Some properties of dis- 
placed squeezed states in the context of second-order coher- 
ence have been investigated before Such states can 
exhibit behavior ranging from photon anti-bunching to super- 
bunching, provided that the state is sufficiently pure, and the 
squeezing weak. Using our modified HBT interferometer, 
and exploiting the high stability and low optical loss of our 
experimental setup, we were able to prepare and measure 
displaced squeezed states that clearly demonstrated photon 
anti-bunching. In addition, we investigated the immunity of 
second-order coherence measurements to optical attenuation. 

Theory: The second-order temporal coherence function al- 
lows the investigation of photon grouping effects in terms of 
normally ordered intensity-intensity coiTelations. The nor- 
malized temporal coherence function of a stationary radiation 
field consisting of a single mode a is defined as the probability 
of a joint detection of two photons at delayed times t and t + r 
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FIG. 2: (Color online) The g*^^'(0) parameter space for a displaced 
squeezed state with squeezing r and displacement a. Regions ex- 
hibiting photon anti-bunching are marked progressively darker, with 
contours giving the precise values. 
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For classical fields the operators are replaced by complex 
numbers, and reordered such that at no delay we can write 
5(2) (0) = with / the field intensity. This quantity 

is bounded below by one, so classical fields cannot display 
anti-bunching. The non-commutation of a and a) prevents 
the quantum expression from being reordered in this way, and 
allows for values of g^"^^ less than one. It is possible to rewrite 
Eq. Q as a product of intensities after a beamsplitter, such as 
measured on the HBT interferometer 
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If 6 = (a + v) /y/2 and c — {d — v)/V2 with V a vacuum 
mode, then Eq. Q and Eq. (|2} are formally equivalent. The 
temporal coherence function witnesses the anti-bunching ef- 
fect when 5 f^^) (q) < g(^)(^T^O). The opposite inequality holds 
true for photon bunching. This method of measuring g'^^) (^-j-'j 
is possible only when one has access to the photon number 
coiTelations directly. 

For experiments in the CV regime, a connection must be 
found for measurements made via homodyne detection. A 
pair of homodyne detectors will measure correlations between 
the field amplitude X+ — a + or phase — — i(d — a^), 
quadratures. Therefore we re-express the coherence function 
Eq. (|3 in terms of quadrature operators 



(3) 



with i, j = +, — and k = b,c. What makes this measurement 
technique possible is that in the equation there are no cross- 
quadrature correlation terms for a single mode. Hence each 
coiTelation or variance term can be measured independently 
by recording the output of the homodyne detectors (H1/H2 in 
Fig. 1), and a value for g*^^) constructed accordingly. 

We wish to consider the coherence function Eq. for a 
weakly squeezed and weakly displaced vacuum state. We 
work in the Heisenberg picture. The operator a is first trans- 
formed according to 

{a)S^ {r)aS{r)D{a) = a + acoshr - sinhr (4) 

where D and S are the unitary displacement and squeezing 
operators respectively fT], and we choose a eSR and r e5R. 
Then by using Eq. Q in Eq. (|3 with t = we find 



5^(0) = 1 + 



sinh^(r) (2q;^ + cosh(2r) - 2a'^ coth(r)) 
(a^ + sinh^(r)) ^ 



(5) 

This function is plotted in Fig. 2 where bunching or anti- 
bunching statistics can be found with the coiTect choice of 
r and a. The stronger cases occur for states approaching the 
vacuum state a = r = for which g'^) (0) is undefined. Indeed 
it's possible to approach the vacuum state while holding any 
value of g^'^\0) constant. The expression for g'^^^O) can be 
extended to arbitrary r by repeating the derivation but includ- 
ing a delay of t on mode b we obtain 
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+ (2 + [a(0) , (r)] + [a(T) , {0)])a\mh'' (r) 

+ (l + [a(0),at(r)]2)sinh*(r)}. (6) 

Here, the commutation relations depend on the shape of the 
filter used in the measurement process. The filter selects a 
particular frequency mode given by 



aiT)=N' 



(7) 



where /(w) is the filter function and N—J^^ /{cjyduj is a 

normalization factor. Using [d^, al,,] ~ 6{uj — lj') and choos- 
ing the filter ^ 57)=1 and zero elsewhere gives the com- 
mutation relations [a(0),a^(T)] = [a(r), a'l'(O)] = sinc(17r). 
Substituting these into Eq. (|6j gives 
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This expression is specific to pure states, but can be gener- 
alized for any Gaussian state having variances V-^^ , V-^ and 
mean ain as measured at the input to the HBT interferometer 
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(see Fig. 1). The general form is 
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This equation was used for making theoretical predictions that 
were based on measurements {V^^; V^^; ain}- A comparison 
with the actual measured g^^^ (r) could then be made. 

Experiment: The schematic shown in Fig. 1 includes: the 
source of squeezed light and displacement operation; the mod- 
ified HBT interferometer with homodyne detection and data 
acquisition. A Nd:YAG laser emitted 1.5 W of coherent light 
at 1064 nm with a line-width of 1 kHz. A major fraction of 
this light was frequency-doubled and used to pump an opti- 
cal parametric amplifier (OPA) composed of a MgO : LiNbOa 
crystal. The OPA could be controlled so as to de-amplify the 
transmitted seed beam by gain factors from 0.9 to 0.5. This 
produced a bright beam (w 5 /iW) having broadband squeez- 
ing of the amplitude quadrature covering 3 MHz up to the cav- 
ity line-width of 15 MHz. Corresponding levels of squeezing 
were measured (at 6MHz) from = 0.89 to = 0.55, 
with values of purity (VJjJ^V^j^) ranging from 1.005 to 1.18, 
respectively. The displacement operation was performed by 
injecting a bright amplitude-modulated field at the 98:2 beam- 
splitter. The relative phase at the beam-splitter was controlled, 
so that the quadrature of the squeezing was preserved. 

With the desired displaced squeezed state prepared, it was 
then subjected to the second-order coherence measurement 
process as performed via the modified HBT interferometer. 
The prepared state was mixed with a vacuum state on a 50:50 
beam-splitter Each of the two resulting output beams were 
sent on toward a balanced homodyne detector, at whose core 
was a pair of matched photodiodes (ETX500) with quantum 
efficiency estimated at 93 ±3%. Care was taken to keep the 
optical path lengths equal. The local-oscillators (LO) were 
mode-matched into the test beam with 96% visibility, and the 
homodyne condition fulfilled by a signal/LO power ratio of 
1/1000. The total detection efficiency was 86%. The LO 
phase was controlled, so as to measure the required quadra- 
ture of either amplitude or phase. The electronic signals of 
the photodetectors were amplified, mixed-down at 6 MHz, 
and low-pass filtered, before being over-sampled by a 12-bit 
digital-to-analogue converter at the rate of 240 kS/s. The time 
series data was then low-passed with a digital top-hat filter 
and down-sampled to 120 kS/s, thereby ensuring a flat power- 
spectrum with bandwidth 60 kHz. (This filter shape deter- 
mines the final g^^Hr) function). The variances and corre- 
lation coefficients of the four permutations of the quadrature 
measurements were calculated using approximately one mil- 
lion data points acquired over 10 successive runs. These val- 
ues were entered into Eq. (|3} to get g'-^^ (t). The uncertainty in 
each measurement (68% confidence interval) was calculated 
using standard methods of error-analysis. 

In this experiment we had the ability to freely vary the 
amount of displacement, and also add broad-band noise to the 
amplitude quadrature which simulated a thermal state. The 
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FIG. 3: (Color online) Experimental measurement of g'^^^r) with 
time delay r in units of bandwidth {rfl). (i) displaced squeezed state, 
(ii) coherent state, (iii) biased thermal state, curves are theoretical 
predictions. 




FIG. 4: (Color online) Main figure and inset: Experimental mea- 
surement of g^'^\0) as a function of displacement ain. (i) displaced 
squeezed state, (ii) weak displaced squeezed state, (iii) coherent 
state, (iv) biased thermal state, curves are theoretical predictions. 



squeezing strength and purity were linked together by the op- 
tical loss and seed-coupling mechanisms of the OPA, but some 
selection was possible by independent variation of the pump 
and seed powers of the OPA. Additional optical loss could be 
introduced by way of a variable optical-attenuator placed be- 
fore the HBT interferometer Finally, the electronic signals 
of each homodyne detector could be temporally shifted with 
respect to one another, both before and after digital sampling. 

Results: Photon anti-bunching statistics from a displaced 
squeezed state were confirmed by the experimental results. 
Fig. 3(i) shows how the measured value of g'^' (0) varied as a 
function of time delay r. In the experiment, a squeezed state 
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was prepared and measured to have {V^'^ = 0.902(1); V^^ = 
1.137(1)}. The state was then displaced by = 0.257(1) 
which was the amount predicted to minimize g^"^^ for those 
particular variances. This prepared state yielded g^'^\0) = 
0.44(22) for zero time delay. Then as r was increased, g^"^^ (r) 
approached unity. Thereby fulfilling the requirement for pho- 
ton anti-bunching i^^^) (0) < ^^^^ (r). 

Displaced squeezed states present their most interesting be- 
havior when (0) is plotted along a range of displacements. 
Here, the state can easily be manipulated into exhibiting either 
anti-bunching or super-bunching statistics, depending only on 
the amount of displacement applied. The experimental data 
is plotted in Fig. 4(i), together with the theoretical curve as 
determined solely by {V^n;l^n}- Super-bunching statistics 
were measured (see inset of Fig. 4) from the prepared state 
{V.+ = 0.901(3); l^- = 1.136(1)} with zero displacement 
= 0.001(2) which produced g'-^\0) = 28(10). On in- 
creasing ttin Still further, the state first found a minimum value 
corresponding to anti-bunching g^^^(O) = 0.41(12) and then 
monotonically approached one. 

Aiming to observe even stronger anti-bunching, we pre- 
pared a squeezed state approaching higher purity {^^n = 
0.890(2) X 1.129(2) = 1.005(3)}. This near pure state was 
produced when the seed power to the OPA was reduced, which 
further de-coupled extraneous noise sources. The g^^^O) 
of this state was measured for a range of displacements; 
Fig. 4(ii). The minimum value was found at = 0.252(2) 
which yielded g^^-'(O) = 0.11(18). For comparison, a two- 
photon Fock state would be limited to g^^^ (0)=0.5. 

To fully characterize our HBT interferometer, we used two 
test states that could readily be produced, and whose prop- 
erties were well known; the coherent state and the biased 
thermal-state (both with variable displacement). Theory pre- 
dicts that a coherent state will produce g'^'^^r) = 1 indepen- 
dent of ain. This was confirmed for r = in the measurement 
Fig. 4(iii), and also for variable r in Fig. 3(ii). Both sets of 
measurements kept within g'^^^(T) = 1.00(6). 

A biased thermal-state was made by taking a coherent state, 
and applying to only the amplitude quadrature, a displace- 
ment driven by a broad-band white-noise source. The state 



was called biased because {Vj^>l; ^^ = 1}. The prediction 
for a strongly biased thermal-state is g^"^^ (0) ~ 3 for V^^ 1 
and small displacments (ain~l). We prepared a state with 
{T^+= 12.80(9); 1^,7=1.039(1)} and varied from zero to 
0.65(1) to obtain Fig. 4(iv). The result was g'-^^O) =2.98(1) 
which adhered to the theoretical prediction. This state was 
also studied under variable time delay; Fig. 3(iii). At t=0, 
the function was maximum g^^-'(O) = 2.98(1), and then fell 
towards one with increasing r, according to the sine function. 
The prepared state was {^^+=14.60(2); yj7=1.025(8); 
0.258(1)}. The results from both the coherent state and 
biased-thermal state measurements agreed with the theoreti- 
cal predictions, and validated our HBT interferometer. 

A curious property of the g'-^^ measurement is that its value 
is unchanged by prior mixing of the state with a vacuum state, 
hence g^"^^ measurements are invariant to optical loss. We 
tested this property by preparing a displaced-squeezed state 
{T^+= 0.894(2); T7,;= 1.139(2); a'i„= 0.255(2)} and pass- 
ing it through variable optical attenuation. With zero atten- 
uation, and total detection efficiency ?7det=86% we measured 
.g'^^^(O) = 0.67(16). Increasing the attenuation to ?]dct=43% 
yielded (?*^^)(0) = 0.43(36). This showed that within a confi- 
dence interval, g^^^ was invariant to optical loss. The invari- 
ance comes at the cost of increasing the uncertainty in the 
measurement. In principle, the total measurement time can be 
increased to accommodate an arbitrary level of attenuation. 

In summary, we have used a continuous-variable measure- 
ment scheme to experimentally probe the second-order coher- 
ence function of quantum states of Ught. The scheme was 
based on the Hanbury-Brown Twiss intensity interferometer, 
but used homodyne detection of the quadrature amplitudes, 
in contrast to discrete-variable single photon counting. By 
preparing displaced squeezed states with the appropriate pa- 
rameters, we were able to observe photon anti-bunching as 
witnessed by 9'^^^(0) = 0. 11(18), which is a non-classical ef- 
fect, and a direct manifestation of the quantum nature of light. 

This work was supported by the Australian Research 
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No. 1 P03B 137 30. 



[1] R. Hanbury Brown, and R. Q. Twiss, Nature 177, 27, (1956). 
[2] M. Henny, S. Oberholzer, C. Strunk, T. Heinzl, K. Ensslin, 

M. Holland, and C. Schonenberger, Science 284, 296 (1999); 

W. D. Oliver, J. Kim, and R. C. Liu, Y. Yamamoto, Science 

284, 299 (1999). 
[3] M. Yasuda, and F. Shimizu, Phys. Rev. Lett. 77, 3090 (1996). 
[4] D. F. Walls and G. J. Milbum, Quantum Optics (Cambridge 

University Press, 1994). 
[5] R. Hanbury Brown, and R. Q. Twiss, Nature 178, 1046 (1956). 
[6] G. Goldhaber, S. Goldhaber, W. Lee, and A. Pais, Phys. Rev. 

120, 300 (1960); D. H. Boal, C.-K. Gelbke, and B. K. Jennings, 

Rev. Mod. Phys. 62, 553 (1990). 
[7] R. Glauber, Phys. Rev. 130, 2529 (1963). 
[8] H. J. Kimble, M. Dagenais, and L. Mandel, Phys. Rev. Lett. 39, 

691 (1977); F. Diedrich and H. Walther, Phys. Rev. Lett. 58, 



203 (1987); G. Leuchs, and U. L. Andersen, Laser Physics 15, 
129 (2005). 

[9] J. G. Rarity, P. R. Tapster, and E. Jakeman, Opt. Commun. 62, 
201 (1987); W. A. X Nogueira, S. R Walbom, S. Padua, and 
C. H. Monken, Phys. Rev. Lett. 86, 4009 (2001). 

[10] M. Koashi, K. Kono, T. Hirano, and M. Matsuoka, Phys. Rev. 
Lett. 71, 1164(1993). 

[11] P. Michler, A. Imamoglu, M. D. Mason, P. J. Carson, 
G. E Strouse, and S. K. Buratto, Nature 406, 968 (2000); 
C. Santori, D. Eattal, J. Vuckovic, G. S. Solomon, Y. Ya- 
mamoto, Nature 419, 594 (2002). 

[12] B. Lounis, and W. E. Moemer, Nature 407, 491 (2000); B. Dar- 
quie, M. P. A. Jones, J. Dingjan, J. Beugnon, S. Bergamini, 
Y. Sortais, G. Messin, A. Browaeys, and P. Grangier, Science 
309, 454 (2005). 



5 



[13] M. Schellekens, R. Hoppeler, A. Pcrrin, J. V. Gomes, D. Boiron, 
A. Aspect, and C. I. Westbrook, Science 310, 648 (2005). 

[14] A. Ottl, S. Ritter, M. Kohl, and T. Esslinger, Phys. Rev. Lett. 
95, 090404 (2005). 

[15] D. F. McAlister, and M. G. Raymer, Phys. Rev. A 55, R1609 



(1997); J. G. Webb, T. C. Ralph, and E. H. Huntington, Phys. 
Rev. A 73, 033808 (2006). 
[16] D. Stoler, Phys. Rev Lett. 33, 1397 (1974); M. H. Mahran, and 
M. V. Satyanarayana, Phys. Rev A 34, 640 (1986). 



